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The non-linear conductance of semiconductor heterostructures and single molecule devices exhibit- 
ing Kondo physics has recently attracted attention. We address the observed sample-dependence 
across various systems by considering additional electronic contributions present in the effective 
low-energy model underlying these experiments. To this end we develop a novel version of the 
superperturbation theory in terms of dual fermions on the Keldysh contour. We analyze the role 
of particle hole asymmetry on the transport coefficients. Our approach systematically extends the 
work of Yamada and Yosida and others to the particle-hole asymmetric Anderson model. It cor- 
rectly describes the strong coupling physics and is free of internal inconsistencies that would lead 
to a breakdown of current conservation. 
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The influence of electron-electron correlation and 
electron-lattice interaction on the flow of charge and en- 
ergy in matter has been a central topic of both explo- 
rative and applied materials science research [H-li. While 
the calculation of transport coefHcients on the basis of 
the fluctuation-dissipation theorem is fairly well devel- 
oped, no generally valid methods exist to go beyond 
the linear-response regime, as e.g. a Boltzmann equa- 
tion based approach relies on well-defined quasi-particles 
and relaxation-time-like approximations. Of particular 
current interest is therefore the effect of strong electron- 
electron correlation on electrical and thermal conductivi- 
ties beyond the linear response regime. Kondo-corrclated 
quantum dots have served as ideal model systems to ad- 
dress this interplay between out-of-equilibrium dynamics 
and strong correlations both experimentally and theo- 
retically. In equilibrium, the Kondo effect leads to an 
enhancement of the linear conductance G = dl /dV\v=o 
to close to twice the quantum of conductance at suffi- 
ciently low temperatures and symmetric coupling (/ is 
the current through the quantum dot and V the applied 
bias voltage) independent of details of e.g. the density 
of states of the leads 0-01 • The fate of this universal- 
ity away from equilibrium has been a subject of intense 
research [sl-ll^l . Recently, the universal aspects of steady- 
state charge transport in the Kondo regime beyond lin- 
ear response through semiconductor heterostructures and 
various single molecule devices have been addressed ex- 
perimentally It was found that the prefactors a 
and 7 of the non-linear conductance, defined via (fcs = 1) 
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differ significantly across different classes of devices. 
Here, T is temperature. Go = G(T 0,0) and Tk is a 



dynamically generated low energy scale, the Kondo tem- 
perature. 

Our primary motivation is to address the discrepancy 
between the results reported in Q {chg = O.I, 7g = 0-5) 



and [15| (as = 0.05, 7s = 0.1) within the single-level 



Anderson impurity model (SIAM) as the effective low- 
energy model for these devices. In the strong coupling 
regime, this model is equivalent to the Kondo model plus 
a potential scattering term generated away from particle- 
hole (p-h) symmetry. As realistic devices are generically 
not p-h symmetric it is important to understand the ef- 
fect of p-h asymmetry on transport properties. 

Although, Kondo physics has been at the center and 
forefront of condensed matter physics not much is known 
about departures from the linear response behavior, as 
coded to lowest order in the coefficients a and 7. A 
full solution of the SIAM out of equilibrium is not avail- 
able and the calculation of these transport coefficients is 
challenging. Results for a obtained from exactly solv- 
able cases are not directly applicable [1, [13] ■ Stan- 
dard approaches, e.g. the numerical renormalization 
group (NRG) yield only linear response transport coef- 
ficients [1^. Selfconsistent methods can in principle be 
extended to the non-linear response regime. They are 
conserving by construction [19| but either fail to cap- 
ture the correct ground state as as in the case of the non 
crossing approximation or the extension onto the Keldysh 
contour is too involved |2d|. Bare perturbation theory is 
in principle applicable [2l| but its non-equilibrium ex- 
tension suffers from an internal inconsistency that leads 
to non-conservation of the charge current away from p-h 
symmetry p^ . The theoretical treatment of the non- 
linear transport within such frameworks therefore as- 
sumes p-h symmetry to avoid issues with current conser- 
vation. In the strong-coupling limit and at p-h symmetry 
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a ~ 0.15 has been obtained independently of the amount 
of asymmetry in the lead-dot coupling between the two 
leads HSliili^. 

As the potential scattering term is a marginally irrele- 
vant perturbation it is expected to modify the transport 
coefficients but its effect should be pertubatively acces- 
sible. To address its effect in a current conserving ap- 
proximation we develop a novel 'superperturbation the- 
ory' on the Keldysh contour for the SI AM. The term su- 
perperturbation theory was first used by Hafermann et 
al. in their treatment of a quantum impurity [27 , 
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small number of bath states defines a reference system to 
set up a superpcrturbative expansion around the numer- 
ical solution in terms of dual fermions. This numerical 
approach has also been extended to the Keldysh con- 
tour [2^. Our approach is rather different: We organize 
a superperturbation theory in terms of the p-h asymme- 
try around the SIAM. This systematically extends the 
work of Yamada and Yosida and others to the asymmet- 
ric SIAM [iil,!!^-!!!. The reference system is solved in 



terms of the renormalized perturbation theory of Hewson 
for the strong-coupling limit of the p-h symmetric SIAM 
extended to the Keldysh contour [iflllll, 34 1- This allows 
us to calculate the transport coefficients analytically and 
give a full characterization of the devices in [3, [lH in 
terms of the SIAM. 

The SIAM Hamihonian is = Hc + Hd + Hd-c, where 
we defined 
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Here, He is the Hamiltonian for electrons in a single con- 
duction band at the metallic leads. Hd is the Hamiltonian 
for localized states in the dot, including the Coulomb in- 
teraction, and Hd-c is the coupling term between the dot 
and the leads. We have defined Ed = ed + U/2. For the 
p-h symmetric case ed = —U/2 and hence Ed = 0. The 
generating functional on the Keldysh contour is given by 



(3) 



where the action on the Keldysh contour is expressed in 
terms of a functional integral over time-dependent Grass- 
mann fields, ipk\a{t) = (cfcAa(0: CfcA<T(0)^ and = 
(dgT (t), Here, the indexes ± refer to the time- 

ordered (-) and anti-time-ordered path along the 
closed Keldysh contour [35|- The lead electrons are non- 
interacting and the resulting Gaussian integrals in Eq. ([3]) 
can be carried out. This results in 
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The effective action S is defined as 
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2tt 



SIV-^ f|>t $ 1 



(5) 



(6) 



is the effective action for a p-h symmetric {Ed 
interacting ([/ 7^ 0) system, and 
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for A = L,R. 
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For simplicity, we assume that the density of states of 
left and right lead, px{ui) — 5{uj ~ e-kx)i are identical 
and p-h symmetric. P-h asymmetry in p\{oj) will have a 
similar effect as the local p-h asymmetry generated by a 
finite Ed- Each lead (L/R) is in equilibrium and is char- 
acterized by temperature (T) and its chemical potential 

(pL/fJ-B)- 



We use the superperturbation scheme [27j to treat the 
term proportional to Ed in the effective action, Eq.(I5|), 
by using as a reference system the p-h symmetric and 
interacting system described by the effective action 
5c/ [$L,^^^] in Eq.®. This yields 

G.,,. = -£rf"'<T3 + (ga^uEda^y^ Gi^^ iEda3g,,u.y' ,(8) 

where G^^^ is the Green's function matrix for the inter- 
acting {U ^ 0) asymmetric {Ed ^ 0) SIAM, 0-3 is the 
third Pauli matrix, and g^^^^ is the Green's function for 
the interacting {U ^ 0) and symmetric (Ed = 0) SIAM. 
Finally, is the dual fermion matrix Green's function, 
obtained from the solution of the matrix Dyson equation 
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and the bare dual fermion Green's function is defined by 
the expression 



EI'^3) 



(10) 



By construction, the expansion works for small and large 
Ed 27|. The explicit expression for the retarded self- 
energy obtained from our superperturbation scheme, up 
to 0{T'^V^), is 



= (1 - x++y + Ed- X+lEd { ^ 
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with X++ = 1 + (3 - 7rV4)(f//^A)2 + 0([/4) [H Hil 
and C = 3/3/(1 + /3)^ where f3 ~ Tl/^r measures the 
asymmetry in the lead-to-dot couplings. It is easy to 
see that our results reproduce the resonant level case for 
{U = 0) and also the special case Ed = U/2 at V = 
considered in [2l| thereby establishing the power of this 
approach. 

We now turn to the evaluation of the transport coeffi- 
cients. The steady-state current through the dot follows 
from the continuity equation [36| 
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where A{uj, V) is the local spectral function (in the pres- 
ence of the dot-lead coupling); /l/ Ir is Fermi function 
in the left/right lead, respectively and V = {fiL — l-J'R)/e. 

We now demonstrate that this scheme gives the cor- 
rect conductance even away from p-h symmetry. Cur- 
rent conservation holds if (i^)S^ (w) — ('^)S^ (lu) ~ 



[22|. This in turn is ensured if Fu{uj) — Fij{u)) 



where the distribution functions are defined by {lu) = 
-F[/(w)(G''(cj) - G"(w)) and S<(a;) = -F[/(w)(E''(w) - 
E°(a;)). It is easy to see that any self-consistent the- 
ory will fulfill this condition [l^. To demonstrate that 
our scheme is indeed current conserving up to the order 
of approximation 0{U^ , E^), we note that the retarded 
component of the self-energy is = — (S^j" + S^~) = 



5^;,. + Ed 



S^'++ so that 
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which implies that this is indeed a current-conserving 
approximation. In Eg. ljlSp we have chosen the chemical 
potentials as to satisfy Tl\ij.l\ = Fijl/ifll for convenience. 
For the special case of identical leads, the distribution 
function of the local Green function is given by /eff (<^) ~ 



{fL{io)rL + fR{Lo)TR)/{TL+rR) 137|. To the same order 
of approximation as the self-energies are calculated, the 
effective distribution function for the interacting system 
is Fuiiu) ^ Uff{u) - (C/12)(C//^A)2(7rreF/A)2/;'(w), 
with fo{uj) the Fermi-Dirac distribution at fi = 0. 
This again demonstrates the current-conserving nature 
of our scheme. The non-linear conductance follows from 
Eq. ((T2|) and the approximation for A{uj,V), obtained 
from the retarded self-energy Eq. pTj) : 
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where 



G{T,V = 0) = Go 
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The renormalized parameters are defined as ed = 
Ed/ A, A = x;+A, u = xll{U/nA), and Gq = 
(2e^/ft-)^(l + (1 — w)^?^)^^ which reproduces the ex- 
act result from Friedel's sum rule up to order O(Me^) as 
[sin(^ — tan~^(ed))] ~ 1— e^. The transport coefficients 
cy, Ctv are to O(e^) given by 



C+y(l + 2C) + e1[5(2^i-l) + (5. 
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CTV = Y [4 + + 2u^ - (2 + - ^4)^ _ ~2 
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In Fig. [U we show our results for a and 7 for various 
cuts through parameter space (ii, id, P). Note, that in the 
strong coupling limit ({t 1) the dependence on id van- 
ishes reflecting the fact that this limit is p-h symmetric 
(sec Fig.[TJa) and (b)). 7 retains its dependence on (3 in 
this limit while a becomes independent of /3 for w — > 1. 
Fig. [TJc) and (d) show the ratio j/a. Eq. ((T7)) shows 
that CvEa and c^ve^ are proportional to the product of 
lead-dot asymmetry /3 and p-h asymmetry id- This prod- 
uct and hence Cve^ and Ctve^ are expected to be small 
in most experimental realizations. Interestingly, the ra- 
tio of CvE^ and Ctve^ is only a function of u. For the 
p-h symmetric case our expressions reduce to the results 
obtained by Oguri and others [k" 



23|,l25|. 



We are now in a position to address the experimental 
results for a = cy /ct and 7 = ctv I Ct mentioned in the 



introduction IJ, [15| . A major experimental challenge is 
to reliably extract the dynamically generated low-energy 
scale A — Tk- The phcnomcnological formula 
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is commonly employed to extract Tk [6| . The parameter 
s apparently fixes ct from Eqs. lfTS]) and (flSl) (s = 0.21 
as in [1^ leads to ct ~ 5.5 and s = 0.22 [l5j results in 
Ct ~ 4.9). It is clear from Eq. that ct is not only 
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FIG. 1: (Color online) (a) a and (b) 7 versus u for various 
id and /3. (c) ratio 7/a versus id and (d) the asymmetry in 
the lead-dot couplings /3 for selected u. In (a)-(d), we used 
the full expressions in id whereas Eqs. (|16fl and (|17p include 
terms only up to order O(e^). 



27| . Our approach is current conserving. This scheme 
may hence be used to characterize nano-structured de- 
vices in terms of the underlying model beyond empirical 
formulas. It can also be applied to obtain magneto- and 
thermal transport properties. Finally, it is possible to 
generalize the results for the single-level Anderson model 
to more complex systems. 
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Note added After completion of this work we became 
aware of Ref. [s^ , which addresses the effect of p-h asym- 
metry on a within a perturbation theory around the p-h 
asymmetric case. A problem with this approach is that 
it fails to recover p-h symmetry at u = 1 and in addition 
it gives an unphysical linear in T term in the spectral 
density away from half filling 71 = 1 . 



a function of u but also depends on the p-h asymmetry 
through id- This complicates the experimental extraction 
of Tk- On the other hand, a and 7 are ratios of trans- 
port coefficients and insensitive to the precise definition 
of Tk- Fig. [TJc) shows the ratio 7/a versus p-h asym- 
metry. From the reported values [IJ, |l5[ it then becomes 
clear that charge fiuctuations (i.e. no p-h symmetry) are 



According to [IJ, [15| the 



present in both experiments. 

coefficients Cve^ and Ctve^ are indeed vanishingly small. 
Yet, they may have been detected in 
mental values reported in 



15|. 



The experi- 
14| are compatible with e.g. 
u = 0.45, ed = 0.1, 13^1 yielding a = 0.1 and 7 = 0.47. 

While we can reproduce as or 75 of [l^ independently, 
it is not possible to reproduce both consistently within 
the SIAM. The underlying low-energy model of the ex- 
periment [iBl cannot simply be the SIAM. One possible 
generalization is that more than one level participates in 
the low-energy properties. Then, already Go is no longer 
given solely in terms of the occupation Ud and the lead- 
to-dot couplings will enter explicitly [ssj . A more likely 
alternative is that local phonon modes renormalize the 
transport coefficients a and 7 differently. 

Our scheme thus constitutes a convenient analytic way 
of characterizing an experimental device in terms of the 
renormalized parameters u, Cd and (3 and the low-energy 
scale A. Inverting the definitions for the renormalized 
parameters will give access to the bare ones 34 1. 

In summary, we have developed a novel analytic 
scheme based on dual fermions to obtain non-linear trans- 
port coefficients for the strong coupling limes of the An- 
derson model. This approach is reliable and controlled 
as demonstrated. It is rooted in the marginal irrelevance 
of the potential scattering term generated by p-h asym- 
metry and on the approaches of 



32l . |34| | as a refer- 



ence system as well as on the dual fermion approach of 
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